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Abstract. Material design at submicron scales would be profoundly affected if the 
formation of precipitation patterns could be easily controlled. It would allow the 
direct building of bulk structures, in contrast to traditional techniques which consist 
of removing material in order to create patterns. Here, we discuss an extension of 
our recent proposal of using electrical currents to control precipitation bands which 
emerge in the wake of reaction fronts in A+ + — > C reaction-diffusion processes. 
Our main result, based on simulating the reaction-diffusion-precipitation equations, is 
that the dynamics of the charged agents can be guided by an appropriately designed 
time-dependent electric current so that, in addition to the control of the band spacing, 
the width of the precipitation bands can also be tuned. This makes straightforward 
the encoding of information into precipitation patterns and, as an amusing example, 
we demonstrate the feasibility by showing how to encode a musical rhythm. 

Keywords: Pattern formation (Theory), Chemical kinetics. Coarsening processes, 
Nonlinear dynamics. Kinetic growth processes (Theory) 
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1. Introduction 

Information encoding and material design involves the creation of patterns which, in 
practice, often means that structures must be produced in a homogeneous media. At 
submicron range which is the target for downscahng of electronic devices, the control 
over the desired patterns becomes difficult and, furthermore, the expenses of traditional 
top-down methods (such as e.g. lithography where material is removed in order to create 
structures) grow steeply. A possible way out of the difficulties is the so called bottom-up 
design where one aims at forming structures directly in the bulk. Nature, of course, 
provides us with illuminating examples of three-dimensional pattern formation at all 
scales [HIS]- Among them, there is a much studied class of reaction-diffusion processes 
yielding precipitation patterns [3], and these processes - suitably planned and controlled 
- are promising candidates for bottom-up designs [H [5] . Indeed, there has been a series 
of attempts to control the emerging patterns by appropriately chosen geometry |6j and 
boundary conditions [5j , or by a combined tuning of the initial and boundary conditions 
[71 [8]. Unfortunately, the above methods of control are not practical enough, and more 
flexible approaches are required. 

Recently, we introduced a novel method of pattern control [9] based on the use of 
electric currents for regulating the dynamics of the reaction zones. We showed both 
theoretically and experimentally that, by controlling the reaction zones, the positions of 
precipitation bands were predesignable. Here we further develop the theory of the new 
method and show that, in addition to the control of the spacings of the precipitation 
bands, it is possible to control the widths of the bands, as well. Thus an extra degree of 
freedom appears which can be used to encode information. We demonstrate the utility 
of this extra freedom by encoding rhythmic patterns into precipitation structures. 

In order to describe the new features of the control by electric currents, we begin by 
a brief description of the properties of the all important reaction zones which provide the 
main input to the precipitation process in the wake of the zone (Secj2]). Next, the effect 
of time dependent electric currents is described and the mathematical details needed for 
the simulations of the process are explained (Secj3]). Finally the idea of how to control 
the width of the precipitation band is introduced and examples of information encoded 
into the widths are presented (SecH]). 

2. Understanding natural precipitation patterns 

The basic idea for the pattern control comes from the observation that precipitation 
patterns are often formed in the wake of moving reaction fronts [21 [3] . The motion of 
the front and its reaction dynamics determines where and when the concentration of 
reaction product crosses a threshold thus inducing precipitation. Consequently, and this 
is the essence of our proposal, control over the precipitation pattern can be realized by 
regulating the properties of the reaction fronts. Guiding reaction fronts and tuning the 
reaction rates in them, however, does not appear to be an easy task. In order to explain 
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how it can be done, we turn to the concrete example of Liesegang patterns [21 fTU] . 
They have been studied for more than a century and a wealth of information has been 
collected about the properties of the front dynamics underlying this pattern formation. 



c 
o 



0) 

o 
a 
o 
o 



a 
a 
b 

rate of reaction 





f 



Figure 1. Left panel: Concentration profiles near the reaction zone. Concentrations 
of the reacting ions (a"*", b~), of the counter ions (a~, b^), and of the reaction product 
Co left in the wake of the front are plotted. The reaction zone defined as the region 
where the rate of C production is significantly different from zero moves diffusively 
{xf ^ and remains well localized during the process. Right panel: Moment of the 
creation of a band in the phase separation process - the concentration c reaches the 
threshold value c* at the position of the front, Xf. 



2.1. Properties of a diffusive reaction front 

Liesegang patterns are characteristic examples of precipitation structures formed in 
the wake of moving reaction fronts [llj. The main ingredients are two electrolytes 
A = {A'^,A~) and B = {B^,B^) which react with reaction rate k according to the 

I k 

reaction scheme + B^ — > C. The reaction product C may participate in further 
reactions but, in the simplest case considered here, it just undergoes a phase separation 
process resulting in an insoluble precipitate provided the local concentration is above 
some threshold |T^. In a typical experiment, the electrolytes are initially separated 
with the inner electrolyte B homogeneously dissolved in a gel column while the outer 
electrolyte A is kept in an aqueous solution. At time t = 0, the outer electrolyte brought 
into contact with the end of the gel column and, since the initial concentration, oq, of 
A is chosen to be much higher than that of B (typically ao/bo ~ 100), A invades the 
gel and a reaction front emerges which advances along the column. The motion of this 
front and the amount of reaction product C left behind the front are clearly important 
factors since they determine the input for the precipitation processes. 

The main features of the reaction front (see left panel of Fig. [1]) are well known 
[T3l and can be summarized in the following three points: 

(i) The front moves diffusively. Its position x/ is given by Xf{t) = ^2Dft where the 
diffusion coefficient Df is determined by the initial concentrations (oq, bo) and by 
the diffusion coefficients of the reagents. 
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(ii) The front is localized. Although the width w of the front is slowly increasing with 
time (w ~ t^^^), it is always much smaller than the diffusive lengthscales (~ t^^^) 
present in the problem. Furthermore, the width becomes negligible for fast reactions 

(W~ 

(iii) The concentration cq of the reaction product, C, left in the wake of the front is 
constant with a value depending on the initial concentrations and on the diffusion 
constants of the reagents. 

Clearly, all these properties are important: (i) tells the location of the front, (ii) 
ensures that the position of the front Xf(t) is precisely given, and (iii) provides the 
amount of C produced at the position specified in (i). 

2.2. Phase separation in the wake of the reaction front 

Once the spatial production of C-s is specified by (i)-(iii), the next step of the 
pattern formation is the phase separation of the C-s. It takes place only if their local 
concentration c is above some precipitation threshold, c > c*. The precipitation pattern 
itself is then the result of a complex interplay between the production of C-s and the 
ensuing phase separation dynamics in the wake of the front. Namely, the experimental 
parameters (ao and bo) are chosen such that cq > c* and, consequently, the front produces 
a precipitation band at the very beginning. This band attracts the newly produced C- 
s from the nearby, diffusively advancing front, thus the concentration c in the front 
decreases below c*. As the front moves far enough, the depletion effect of the band 
diminishes and the c > c* condition is satisfied again in the front, thus leading to the 
formation of a new band. A quasiperiodic reiteration of the above process yields the 
Liesegang patterns (see right panel of Fig. [1]). Depending on the details of the phase 
separation dynamics, the position of the n-th band x„ may vary, but there are three 
well established laws that govern the structure of the Liesegang bands: 

(a) Time law [15]: The position of the n-th band x„ (measured from the initial interface 
of the reagents) is given by Xn = ^j2Dftn, where tn is the time of creation of the 
band. 

(b) Spacing law [16]: The positions of the bands form a geometric series Xn ~ (1 +p)"' 
with a spacing coefficient p > 0, such that distances between successive bands 
increase with the band index n. 

(c) Width law [17j: The width of the n-th band Wn is proportional to its position: 

Wn ~ Xn- 

The above laws can be derived [12] by using the Cahn-Hilliard equation [18] for 
describing of the phase-separation process. This approach also allows to demonstrate 
that the band positions can be controlled by oq or bo since p in the spacing law depends on 
these quantities (p ~ l/ag is the so called Matalon-Pacter law [12]). Unfortunately, the 
possible changes are rather limited since the band positions invariably form a geometric 
series. 
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3.1. Main idea - controlling the motion of the reagents by electric current 

Equipped with the understanding of both the front motion and the precipitation 
processes, we can start to think of possible control mechanisms. There are basically 
two ways to change the structure of the pattern characterized by the spacing law (b) 
and the width law (c). First, one can try to change the functional form of the time 
law (a). This can be done by using various geometries or patterns in the initial state 
O El [71 [20], by employing guiding temperature- or pH fields or by considering 
systems where the diffusion of the reacting species is anomalous [2T],[22]. Unfortunately, 
these methods are rather unwieldy and are not flexible enough to easily create arbitrary 
patterns. 
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Figure 2. Experimental setup for producing Liesegang precipitation patterns as 
described in Sec. 12.11 The controlling agent is the generator providing electric current 
I[t) with a prescribed time-dependence (figure is taken from [5]). 



The second method keeps the time law unchanged and aims at controlling the 
creation time t„ of the n-th band. Recalling that t„ is the instant when the concentration 
of C-s crosses the threshold value c*, one recognizes that tn can be controlled by 
regulating the concentration cq{x) at the front. Recently, it was shown both theoretically 
and experimentally that the above method can be made to work by sending a time- 
dependent electric current through the system [9]. The schematic setup for control is 
shown in Figj2] and, at a phenomenological level, its working can be understood rather 
easily. Indeed, consider an imposed current which drives the reacting ions towards the 
reaction zone (we shall refer to this current as forward current). It is clear that the 
forward current enhances the production of C-s in the reaction zone. Reversing the 
direction of current {backward current), on the other, hand works against the reaction 
and results in a lower production of C-s. Thus, provided the position of the front is 
known [i.e. Xn{t) is available from the time- law], the times t„ of crossing of the threshold 
concentration and, consequently, the positions of the band Xnitn) can be controlled by an 
appropriately chosen current. Since managing the electric current is not an experimental 
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difficulty, the above method provides us a flexible technique for the creation of complex 
precipitation patterns. 

3.2. Mathematical description of the process 

For a more quantitative description of the control-by-current method, we shall use a 
mean field model that has been developed in a series of papers during the last decade 
P, [12], [231 [21] ■ The ffist part of the model addresses the irreversible + B~ — > 
C reaction-diffusion process for totally dissociated electrolytes A = {A^, A^) and 
B = {B^,B~) that are initially separated in space. The evolution equations for 
the concentration profile of the ions a^(x, t) and t) are obtained by assuming 
electroneutrality on the relevant time and lengthscales [23] and, for the case of 
monovalent ions with equal diffusion coefficients, the equations are as follows |9| 



Here D is the diffusion coefficient of the ions, j{t) = I{t)/A is the externally 
controlled electric current- density flowing through the tube of cross section A, and 
S = g(a"'" + a~ + b'^ + b~) with q being the unit of charge. The reaction rate k is taken 
to be large resulting in a reaction zone of neghgible width. Note that this assumption is 
compatible with the typical reactions used in experimental setups producing Liesegang 
structures. 

The second part of the model explains the pattern formation through the separation 
of the reaction product C into a high- and low-concentration phases. The evolution of the 
concentration c{x, t) is obtained from the Cahn-Hilliard equation with the addition of a 
source term corresponding to the rate of the production of C-s {ka'^b~) [121 [2l]. The free 
energy driving the phase separation is assumed to have minima at some low (q) and high 
(c/i) concentrations of C and, furthermore, it is assumed to have the Landau- Ginzburg 
form in the shifted and rescaled concentration variable m = (2c — — ci)/{ch — ci). In 
terms of m, the equation describing the phase-separation dynamics takes the form [T2] 



Here S{x,t) = 2ka'^b^ / {ch — ci) is the source term coming from equations ([IHll)- The 
two parameters A and a are fitting parameter at this stage, they can be chosen so as to 
reproduce the correct experimental time- and lengthscales [251 12^- 

Equations ([TM]) are a closed set of equations for the concentrations of the 
electrolytes (a^, b^) and of the reaction product c. Together with the specification 
of the initial- and the boundary conditions, they provide the mathematical formulation 
of the problem. Below we shall consider numerical solutions of these equations which 
were obtained by the classical fourth-order Runge-Kutta method. 



dtb- 
dta' 
dtb+ 



Ddla^ - ]{t) 9^(a+/S) - ka+b 
Ddlb- + ]{t) 9^.(r/S) - ka+b 

Ddla- +j{t)dM'/^) 

Ddlb^-j{t)d,{byi:). 



(1) 

(2) 
(3) 

(4) 



dtm = — AA(m — + a Am) + S{x, t) . 



(5) 



Encoding information into precipitation structures 



7 



3. 3. Simulation results - properties of the front in the presence of a current 

In order to obtain a more detailed understanding of the front dynamics, we studied the 
numerical solution of eqs. ([T]) with initial conditions of separated electrolytes [a^(x < 
0,t = 0) = ao, a±(x > 0,t = 0) = 0, h^{x < 0,t = 0) = 0, h^{x > 0, t = 0) = 6o], and 
monitored both the position of the front Xf{t) and the rate of production S = kab of the 
C-s (a brief account of these simulations has appeared in [9]). The physical parameters in 
the equations were chosen to be close to the experimentally relevant values (ao/&o = 100, 
D = 1.22 ■ 10"^ m^/s, a = 10~^ m^, A = 0.17 ■ 10~^ m^/s) and we considered following 
scenarios for the current. The no current case was used to reproduce the known front 
properties (i-iii). Constant forward and backward currents of amplitude \j\ = 10 A/m^ 
were simulated to check whether the time law holds on the experimentally relevant 
timescales. Finally, we studied the most interesting case of an alternating current of 
constant absolute value (|j| = 10 A/m^) with the sign of it changing in a square wave 
pattern at times with ?t, = 0,1,2,..., and r fixing the timescale of the protocol. 




Figure 3. Left panel: The posftion of the reaction front (measured in units of the 
length L of the gel column) versus dilfusion length in the absence of a current, displayed 
for constant forward or backward current, and for a quasiperiodic current (changed at 
times rn^). Right panel: Concentrations of the reaction product in the wake of the 
front for the cases considered on the left panel [9]. 



The results are shown on FigJSl The left panel displays the front motion and shows 
that in all cases the diffusive nature of the front is hardly changed, i.e. the time law 
(a) remains valid. On the other hand, as can be seen on the right panel of FigOl 
the C-production at the front is strongly influenced by the character of the applied 
current. Backward current leads to a decrease of the concentration of C-s left behind 
the front, and co{x) reaches values below the phase separation threshold co(x) < c* thus 
eliminating the possibility of precipitation. In contrast, the forward current increases 
the C-production steeply and brings the system quickly into the unstable regime thus 
inducing precipitation. It follows then that, in the case of a quasiperiodic current, the 
phase separation can be triggered and timed by switching on the forward field. 
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4. Pattern design 

4.I. Controlling the spacing and the width of the bands 

Once the front dynamics, summarized in Figl3|, is understood, one can invent appropriate 
current dynamics that results e.g. in equidistant band patterns [9]. The wavelength d of 
the periodic pattern can be predesigned by switching on the forward currents at times 
tn = {2n)^T, where n = 0, 1, 2, . . . and r = d'^/SDf. If the desired period d is smaller 
than half the local wavelength of the Liesegang pattern which would be present without 
the current (see Fig{T]), then spurious bands may appear due to a natural increase of 
the concentration of C-s. This can be avoided, however, by switching on the backward 
current when the front is halfway between x„ and Xn+i, i.e., at times {2n + 1)^t. 

One can also create more complex patterns both experimentally and theoretically 
[9]. When creating patterns of several wavelengths, variable widths of the bands may 
also become an important part of the patterns. The issue of the width control, treated 
below, is the novel aspect of the present paper. 

Let us begin by finding an estimate of the width of the equidistant bands having a 



period d = JSDjt. From simulations we know that even in the presence of current, the 



position of the front is well approximated by Xf{t) = ^2Dft where -D/ is given by [13]: 

D, = 2DLr^\^-^^ll^]]\ (6) 

For a typical ratio of the initial concentrations of A and B, ao/bo = 100, this yields 
Df = 5.43-D, where D is the diffusion constant of the ions. 

An important point now is that simulations of the equidistant case [i.e. when the 
forward currents is switched on at times t„ = (2n)^r] reveal (see FigH]) that, although 
the C-production varies strongly within a period, the average concentration is practically 
equal to the zero current case. This means that, for the estimation of the width, we can 
replace the complicated function Co(x) by the result of the homogeneous production of 
C-s 11 



For ao/bo = 100 this yields Co/oo ^ 1.145 x 10^^. 

Using now the conservation of C-s, namely that the amount of C-s produced in a 
period {d cq) is distributed into the high- and low-concentration regions of length w and 
d — w, the width of the bands, w, is calculated as 



w='-^d='-^J^, (8) 

Ch-Ci Ch-Ci^ 

and we find that the width depends on r in the same way as the period. The width also 
depends on the diffusion constant of the front, Df which, in turn, depends on the ratios 
of initial concentrations ao/&o as well as on the diffusion coefficients of the electrolytes. 
Most easily, however, the width can tuned by the initial concentration of A, which is 
proportional to Cq. 
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Figure 4. C-concentration in the wake of the front for the equidistant pattern (red 
dashed hne) and its average over a region of distance d (red soHd Une) in comparison 
with the no-current case (blue line). 



In Figj5] one can see the numerical verification of tlie above approximations for tlie 
width and the wavelength of the periodic pattern. The smaller the value of r the better 
is the estimate, since for large values of r we enter a regime where the deviations from 
the time-law (a) become significant. 




Figure 5. Left panel: Definition of the width w and the period d of the equidistant 
pattern. Right panel: Dependence of the numerically obtained widths (red cross) and 
periods (blue cross) on r in comparison with the theoretical prediction (solid lines). 



4-.2. Encoding information into precipitation structures 

Since spacing and width can be controlled by properly designed current, more complex 
structures should be realizable. The pattern shown in Fig. [6] has been created using 
a slightly generalized method of control described above for the periodic pattern. The 
underlying pattern in Fig. [6] is an equidistant pattern but the width of certain bands 
is increased. This has been achieved by a longer time interval for the forward field for 
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the wide bands. More precisely, in order to make the width of the n-th band larger, 
the switching protocol is changed: instead of switching on the backward current at time 
(2n+l)^r, it is switched on at time ((2n+/)^r such that the band appears approximately 
a factor / larger. Similarly it is also possible to increase the spaces between the bands 
by increasing the time duration of the backward current. In Fig. [7] we demonstrate this 
by creating a structure (the Morse code of the famous number 42 [26]) where a current 
dynamics had to be designed which yields simultaneous control of the spacing and the 
width of the bands. 




Figure 6. Musical rhythm written into Liesegang bands. The upper panel shows the 
concentration-profile of narrow and wide bands (the first four bands are cut out of 
the picture as well as the final ones since they don't belong to the implementation of 
the rhythm). This is compared in the lower panel to the Bolero rhythm composed by 
Ravel ES . 



c(x) 



















j — \ 












' — ' 








f \ 


































u 











X 



Figure 7. Concentration profile, c(x) for a numerical test of the combined width and 
space control. In the simulation the number "42" is written into precipitation bands 

based on the Morse alphabet: 4 ( ) 2 ( ). The dots are implemented by narrow 

bands, the underscore by a wide band approximately a factor f — 2 larger, and the 
characters are separated by spaces, widened by the same factor /. 
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5. Conclusions and outlook 

We have developed further the techniques of the engineering of precipitation structures 
using electric currents. It was shown that an appropriately designed time-dependent 
current allows to control not only the band-spacing but also the width of the bands 
in one-dimensional structures. This permits the encoding of information by controlling 
either only the the width (as exemplified by the Bolero rhythm) or by combining the 
width and space control (as was shown on the example of the famous number 42). 

A naturally arising question concerns the limits of the approach when trying to 
downscale the structures. As far as the timescales are considered, since the diffusion 
coefficients of the participating agents are of order 10~^ m^/s, imposed currents on the 
time scale of r = 0.1s provide control on lengthscales of 10/i. Thus scaling the dynamics 
of the current does not pose an experimental problem. 

The real problem with downscaling is the width of the bands. The minimum width 
is clearly related to the width of the reaction front which is not negligible on the scale of 
microns and, depending on the reagents, can be even at the scale of mm-s [28J . It should 
be noted however that spontaneous pattern formation on the nanoscale level has been 
observed [29], and these type of systems might be appropriate for the pattern control 
by imposed electric currents. Another problem with the width of the bands is that the 
thermal fluctuations in the concentrations of the reagents combined with those of the 
gel make the reaction front uneven and, depending on the surface tension of the created 
bands, they may lead to a roughening of the band on a scale that is comparable to the 
width. Clearly, advances in downscaling can be achieved only if models are developed 
which can treat all the above problems. 

Finally, we note that combining the proposed technique with the already existing 
indirect control strategies such as the choice of geometry and initial conditions opens up a 
wide spectrum of possible structure design. So far the feasibility of the approach has only 
been shown for a few types of predesigned one-dimensional patterns. A straightforward 
extention would be the creation of rings or spheres with a predesigned internal structure. 
Further, the extention of the current-control technique to stamping methods [5j on the 
mesoscopic level could be used to create even more complex pattern designs which may 
become useful in engineering applications. 
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